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Abstract 

This is the second paper in a series on light scattering from optically anisotropic scatterers 
embedded in an isotropic medium. The apparently complex T-matrix theory involving mixing 
of angular momentum components turns out to be an efficient approach to calculating scattering 
in these systems. We present preliminary results of numerical calculations of the scattering by 
spherical droplets in some simple cases. The droplets contain optically anisotropic material with 
local radial or uniform anisotropy. We concentrate on cases in which the scattering is due only 
to the local optical anisotropy within the scatterer. For radial anisotropy we find non-monotonic 
dependence of the scattering cross-section on the degree of anisotropy can occur in a regime for 
which both the Rayleigh and semi-classical theories are inapplicable. For uniform anisotropy the 
cross-section is strongly dependent on the angle between the incident light and the optical axis, 
and for larger droplets this dependence is non-monotonic. 
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I. INTRODUCTION 



In the first paper in this series 0, we developed a T-matrix theory of light scattering by 
spherical but optically anisotropic - either radially or uniformly - scatterers. Whereas for 
radial anisotropic scatterers it is possible to derive closed form expressions for the elements of 
the T-matrix, in the uniformly anisotropic case this is no longer true. To further complicate 
matters, in this latter case the T-matrix elements can only be derived numerically by 
inverting a rather difficult set of matrix equations. Because the spherical symmetry of a 
usual scattering problem is reduced in this case to cylindrical symmetry, the scattering 
involves angular momentum mixing. The consequence is that the set of equations to be 
inverted is in principle infinite in number. 

In this paper we continue this study. We find that notwithstanding the daunting nature 
of the formal structure of the T-matrix theory, in fact this theory provides a viable and 
efficient strategy for calculating the elements of the T-matrix, and thus computing the 
scattering properties of these objects. Throughout this paper we shall refer to the first 
paper in this series as I, and write Eq. (n) of I as Eq. (n.I). 

A number of approaches are available to study light scattering by complex objects. A brief 
summary is as follows. The scattering amplitudes can be described using Green's function 
techniques 0], but these involve solving complex integral equations over infinite domains. 
Under some circumstances one can approximate the kernels of these equations either as the 
incident wave or as a semi-classical perturbed wave, leading to the well-known Rayleigh- 
Gans (R-G) and Anomalous Diffraction Approximations (ADA). These have been used by 
Zumer and coworkers to examine the problems we consider in this paper [Q, ||. However, 
the approximations are only valid over certain wavelength and optical contrast regimes. 
The century-old Mie strategy and its modern T-matrix extensions yield exact solutions, 
but unfortunately this approach does not work in every case. Finally one can of course 
use real space finite element approaches, but these are notoriously inefficient at reproducing 
known solutions. For a more comprehensive review we refer the reader to Chap. 2 in || and 
references therein. 

The T-matrix theory is known to be a computationally efficient approach to study light 
scattering by nonspherical optically isotropic particles [[J. One may thus expect that a 
T-matrix approach to geometrically spherical but optically non-spherical scatterers can at 
the very least enable scattering properties to be evaluated when the approximate methods 
cannot be applied. In addition, whereas the region of validity of the approximate methods 
such as R-G and ADA in the case of isotropic scatterers is reasonably well-understood, in 
the case of anisotropic scatterers this problem has not been been studied in any detail. 

In I we have discussed composite spherical scatterers, consisting of a central isotropic 
core plus a surrounding annular layer in which the optical tensor is anisotropic: = 
e±[Sij + u(h^h)ij], where u = (ey — e±)/e± is the anisotropy parameter. For radial anisotropy 
the optical axis, n, is directed along the radius vector, n = f , and for the uniform anisotropy 
the optical axis is parallel to the z-axis, n = z. These cases present different mathematical 
challenges to the theorist. 

Light scattering from the radially anisotropic annular layer was first studied long ago 
by Roth and Digman || using the technique normally known as Debye potentials. In 
an earlier paper we have recovered this solution as a special case of a more general 
set of anisotropies. A crucial step in the derivation of this result involves writing the so- 
called modified T-matrix ansatz (see Eqs. (26.I)-(27.I)). The spherical symmetry of the 
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scatterer requires the modes (27.1) entering the ansatz (26.1) to be proportional to the 
corresponding vector spherical harmonics. As a result the scattering does not mix different 
angular momenta. The T-matrix is then diagonal over the angular momentum indices and 
the azimuthal numbers and the elements of the T-matrix are expressible in closed form (see 
Sec. IV. A of I). 

The uniform anisotropy case is much more difficult. The light scattering problem for a 
uniformly anisotropic spherical scatterer is not exactly solvable |5| . As an alternative, it has 
been studied by using R-G and ADA g f§. 

In I we approached this problem by formulating it as a suitably modified T-matrix 
theory. It is necessary to relate the plane wave packet representation to expansions of 
electromagnetic fields over vector spherical harmonics. The net result is to define a set of 
wave functions representing exact solutions of Maxwell's equations in an anisotropic layer 
that are at the same time deformed spherical harmonics. The coefficient functions that 
enter the expressions for these wave functions describe angular momentum mixing and are 
computationally accessible. 

In this paper we find numerical results for the total scattering cross-section in the limiting 
case of a droplet, i.e. when the radius of the isotropic core of the scatterer, R 2 , is negligible 
(i?2 0). The scattering by a droplet presents fewer technical difficulties than scattering 
by the annular layer which was our principal focus in I. Anisotropy effects are our primary 
concern and for this reason we pay special attention to the case where the scattering can 
be solely attributed to the presence of the anisotropy. Specifically, in our subsequent cal- 
culations we consider the special case for which the ordinary wave refractive index and the 
refractive index of the material are matched, n = n Q or e = e±. 

The paper is organised as follows. In Sec. II we use the general T-matrix formalism 
developed in I to write down some necessary formulae relevant to the special case of droplets. 
In Sec. [TTT] we make brief comments on the numerical strategy and present some numerical 
results. Finally in Sec. [TV] we make some brief concluding remarks. 

II. T-MATRIX CALCULATIONS FOR DROPLETS 
A. Notation 

In this section we adapt the key theoretical relations derived in I so that they can be 
used in the case of droplets. In addition, we rewrite the expressions for the total scattering 
cross-section in a more convenient form. 

We first introduce some notation. Radially anisotropic droplets present an isotropic face 
to the world, to the extent that the scattering properties are independent of the the direction 
of the incident wave. In the case of the uniformly anisotropic droplets, this is no longer true. 
In this case the scattering geometry is shown in Fig. |1| The angle of incidence 9i nc is the 
angle between the direction of incidence k and the direction of the uniform anisotropy. 
The y direction is perpendicular to the plane made from these two directions. We shall 
show explicitly that the scattering process does not involve the y-component, E y inc \ of the 
incoming plane wave 

Ej nc — ^ exp(z kj nc • r) , ^ = E^ ' 'e^k^c) -|- Ey ^ey(kj nc ) , kj nc = k\s.i nc , 

(1) 
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FIG. 1: Uniformly anisotropic droplet. The polarisation vector e y (ki nc ) is normal to the plane of 
the picture, e y (kj nc ) oc k inc x z. Inside the droplet plane waves linearly polarised along e y (k inc ) 
represent ordinary waves. 

provided the refractive indices n = and n Q = ^/ej are matched. In this paper we put 
the magnetic permittivity equal to the unit and slightly change the notations: n and e± 
instead of n\ and e 1 . This corresponds to the physical condition in which there is a matching 
condition for the refracted ordinary wave inside the scattering droplet. 

B. Equations for T-matrix 

Since the electromagnetic field must be regular at the origin, the harmonics inside the 
droplet are now given by 

Ej m = OLj m Pj m (po, f) — n Q Q.j m Pj m (p , f ) , (2a) 
Hjm = afl Q$™ } (p , r) + n Q afl Q^(p , f ) , (2b) 

where p Q = m kr = m p and m = n /n is the optical contrast. 

The modes and have been given by Eqs. (27.1) for radial anisotropy and by 

Eqs. (43.1) for uniform anisotropy. 

Following Eqs. (32.1), the continuity conditions at the outside of the droplet, r — R±, can 
then be written in matrix notation as follows: 

E( (c) \ / (mc)\ . / n(sca)\ 

a j'm | _ -pj I a jm \ , -pi [ Pjm \ (o\ 



4 



( Up) o 



\ 







-n 



-i 



f Ur) 



[Jj(p)]'J 



I hf{p) 
n[hf{p)\' 


V o 



o \ 



hf\ P ) 
-n-'[hf\p)]'J 



(4) 



where the index 1 indicates that matrix elements are calculated at the boundary of droplet, 
r = R x . 

In the case of a radially anisotropic droplet the matrix on the left hand side of Eq. (|3|) is 
diagonal over angular momentum numbers 



R»-> m {r) = 6 jr 



( jj(Po) \ 

n [jj{po)]' 

]]{po) 

\ -n-%(p )}') 



jO + l)=j(J + !)/(! + «) 



(5) 



When the droplet is uniformly anisotropic, the matrix R> 3 '' ,m {r) is no longer diagonal 
over angular momentum quantum numbers. In this case it takes the form: 



R" 



iPo) 



n q%? m (Pa) 
n q-P; e 2(po) 

\ Pjj';m\Po) 



-1 (m,e) . 



o {e ' m '(o ) 
'iii'-.m \P°) 



-n 



-1 ^( e ' e ) 



(6) 



" Pjj';m(Po)/ 



where the coefficient functions are given by Eqs. (C3.I)-(C10.I) in Appendix C of I. 
The system (B) can be then simplified by multiplying both sides by the matrices 



H j (r) 



-ip 



Uj(p)]' -n 1 j j {p) 



















n' 1 [Up)}' Up) 



£p(r) = -ip 2 



'[hf\p)]' -n-'hf\p) 









,(1), 



■nr 1 [hf{p)} , hf ) {p) f 

Using the fact that the Wronskian for spherical Bessel functions is given by 

W{j J (p) 1 hf\p)} = z/p\ 
we derive a system equivalent to Eqs. (Bt) in the following form: 



(7a) 
(7b) 



(inc) 

a- 

— 1 ~ (inc) 



(sea) 
jm 



n 



(sea) 
jm 



where Ajj>. m 



-U( ■ R{ 



^2 B ^'';' 

j'>\m\ 



j'>\m 



■jj';m 



0£ j'rn / 
a j'm. 



JJ :>n and Bj/:,,, 



ni r{ 



30 ; m 



(8a) 
(8b) 
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These are the equivalent of the equations (51.1) for scattering by a spherical annulus. The 
equations are considerably simpler, as a result of what might be called mode decoupling. 
The crucial point is that the normal modes inside the droplet are all regular at the origin, 
and thus two types of modes which appear in the case of the annulus do not appear here. 
This reduces the number of variables in the problem by one half, even in the case of the 
uniform anisotropy which, at least in principle, presented a certain number of problems in 
the annular case. 

From Eq. (§) we can immediately derive an expression for the T-matrix: 

Tjjv. m = [A ■ B ] 3 y ;m • (9) 

For radial anisotropy all the matrices on the right hand sides of Eqs. (§) are diagonal. 
So, it is easy to write down the result for the T-matrix: 



where 



'T 11 



rpll = [jjiPo)]! \jj(p)]l ~ m Q [jj(po)}'l [Jj(p)]l / Ua) 

' " MPo)\i[hT{p)]'i-mobAPo)\i[hT{ P )]i' 



m Q [;j-(p )]i [jj(p)]'i ~ \j](Po)]'i [jj(p)} 



rpTZ » w "'^ w /JJ. IVJ V V f lib") 

These formulae bear close resemblance to the well known Mie expressions ||. 



C. Scattering efficiency 



We have seen in I that the total scattering cross-section for the uniformly anisotropic 
scatterer depends on the angle of incidence 9i nc . The scattering cross-section also depends 
on the polarisation of the incident wave. In order to emphasise this, let us express the 
coefficients of the incident wave and (see Eq. (9.1)) as follows 



(inc) _ (inc) E (ine) (inc) E (inc) fl^ 

~(inc) _ ~(inc) E (i nc ) --(inc) „( inc ) x^jx 

where 



= V Wi + l)] 1/2 ^(k mc ) , a<£> = [vr(2j + l)] 1 ^^ , 
a&l = [t(2j + l)] 1/2 ^(k mc ) , a<£> = -ni' [7r(2j + l)] 1 / 2 ^^) . (13) 
From Eq. (19.1) the scattering efficiency then can be rewritten in the following form: 
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Angle of incidence 



FIG. 2: Scattering efficiency of uniformly anisotropic droplet as a function the angle of incidence 
(the angle between the incident wave and the optical axis) at various values of the anisotropy 
parameter, u = (eii — ej_)/e_i_, with kR\ = 4.0 and n = n Q . 



s 



n — ^ sca — t- 1 n p(* nc ) Tp( inc )* 

V - ^2 - hnc 2^ 1^ ^aP^a h f3 
1 a=x,y P=x,y 



jm 



Qap-{Kttl) 7T > , [P.7m;aP,m:/3 + n Pjm;aPjm; 



Pjm; a 
-1 3(«») 



(inc) 

OL t 
J m; a 

-33';m J _j ~(inc) 

n a •/ 



(14) 
(15) 

(16) 



Note that for radial anisotropy Q Q /3 = S a pQ- By contrast, in the case of uniformly 
anisotropic droplet we have rather strong dependence of the scattering efficiency on the 
polarisation of the incoming wave. When the refractive indices are matched, n = n D , it is 
expected that the scatterer does not change the y component of the incident wave, which 
simply transforms into the ordinary wave inside the droplet without being affected by the 
scattering process. The algebraic interpretation of this fact is that the amplitudes of the 
scattered wave flj^l and h^fm^y are equal to zero. However, it is not straightforward to 
see that the system (§) is consistent with this conclusion. We show this in Appendix 
the proof involves using relations (46.1) in which the electric and magnetic fields inside the 
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FIG. 3: Dependence of the scattering efficiency on the angle of incidence for uniformly anisotropic 
droplet at various values of the size parameter and u = 0.4. The refractive indices n and n are 
matched. 

droplet are expressed in terms of appropriate normal modes. 
III. NUMERICAL RESULTS 



In this section we present numerical results for the scattering efficiency defined by 



Eqs. (|14]) and ([To]) . We are primarily interested in anisotropy-induced scattering. In or- 
der to concentrate on this test case, we consider the case when the refractive indices n and 
n a are equal. We shall present a more comprehensive analysis of all possible cases, including 
the results for the angular distribution of the scattered waves, elsewhere. We begin with brief 
comments on numerical procedure and then proceed with the description of the calculated 
dependencies. 

It is rather straightforward to perform computations for radially anisotropic droplets. 
The expressions for the elements of T-matrix are known and given by Eqs. ( |lla| )-( pT5| ). We 
can thus evaluate the scattering efficiency by explicitly computing the sum in the expres- 
sion (22.1). 

For uniformly anisotropic droplets T-matrix can only be computed numerically by solving 
the system of equations (|8|) jT(J . Some highlights of the results are presented below. 

The dependence of the scattering efficiency on the angle of incidence is shown in Fig. 0. 
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FIG. 4: Scattering efficiencies of radially and uniformly anisotropic droplets versus the size 
parameter at u = 0.4, 0j nc = tt/2 and n = n Q . 

If the size parameter, kRi, is not very large, the scattering efficiency Q xx is a monotonically 
increasing function of the angle of incidence, 8i nc , in the region from to tt/2. By symmetry 
Qxx(9inc) = Qxxijt /2 — 9 inc ), and so the scattering efficiency decreases in the range from tt/2 

tO 7T. 

In Fig. we show what happens for shorter wavelength and thus higher values of kR\. 
Now, for for relatively large values of the size parameter, the cross-section dependence on 
the angle is incidence is no longer monotonic. For example, at kR\ = 20.0, the angle at 
which the scattering efficiency Q xx reaches its maximum value is no longer at tt/2. 

Fig. [| shows the scattering efficiencies Q xx (tt/2) (for uniform anisotropy) and Q (for radial 
anisotropy) versus the size parameter. The scattering efficiency of uniformly anisotropic 
droplet has a pronounced peak located at about kRi w 10.0 and exhibits strongly non- 
monotonic behaviour. By contrast, the corresponding dependence for the radially anisotropic 
droplet is monotonically increasing. In this case the first maximum is reached at kR\ ~ 20.0, 
outside the range of kRi shown in Fig. [|. 

The scattering efficiencies as a function of the anisotropy parameter, < u < 1, at 
different values of the size parameter are plotted in Fig. ||] and Fig. || for the cases of radial 
and uniform anisotropies respectively. In both cases an increase in the size parameter leads 
to the appearance of peaks in this range of u. As compared to radially anisotropic scatterers, 
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Anisotropy parameter 



FIG. 5: Scattering efficiencies of radially and uniformly anisotropic droplets versus the anisotropy 
parameter at various values of the size parameter for 9i nc = ir/2 and n = n . 

the uniformly anisotropic droplets seem to be more sensitive to changes both in the size and 
in the anisotropy parameters . 



IV. CONCLUSION 

In this work we have described some of the numerical results calculated using the T- 
matrix theory developed in paper I. In particular, we have studied the scattering efficiency 
of radially and uniformly anisotropic droplets in which the ordinary refractive index matches 
the refractive index of the material surrounding them. 

The assumption in which the ordinary refractive index of the droplet matches the isotropic 
dielectric constant in the surrounding medium is not taken in order to simplify the numerical 
treatment. Rather in this paper we wish to study the light scattering properties which can be 
solely attributed to the anisotropic part of the dielectric tensor. Thus we have the anisotropy 
effects separated out to concentrate on differences between isotropic and anisotropic optical 
axis distributions. 

Clearly, the difference in symmetry causes the most crucial differences in the light scat- 
tering. For uniformly anisotropic droplets the scattering efficiency depends on the angle of 
incidence and the polarisation of incoming wave, whereas for radially anisotropic scatterers 
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Anisotropy parameter 

FIG. 6: Dependence of the scattering efficiency on the anisotropy parameter for radially anisotropic 
droplet at various values of the size parameter and n = n . 

it does not. In other words, the scattering from radially anisotropic droplet shares some 
features of isotropic scatterers. 

Nevertheless, the angular dependencies for the scattered wave intensity and the depolar- 
isation factor shown at the end of I clearly indicate the pronounced differences. 

The graphs plotted in Figs. f|-[] indicate that uniformly anisotropic droplets are more sen- 
sitive to changes in the wavelength and anisotropy parameters than are radially anisotropic 
droplets. Our results are also consistent with results of previous studies ||, 0] that the inter- 
nal spatial distribution of the optical axis is a factor which strongly affects light scattering 
from anisotropic scatterers. 

The results of this work can be regarded as the first step towards more comprehensive 
study of light scattering by anisotropic scatterers. We have demonstrated that the T-matrix 
approach developed in I can be used in an efficient numerical treatment of the scattering 
problem. It is thus reasonable to expect that further progress can be made by applying 
this theory to more complex problems involving light scattering by optical anisotropic liquid 
crystalline systems and other related problems, as discussed further in the last section of I. 
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APPENDIX A 



In this appendix we show mathematically, suing our formalism, the physically obvious 
result that if the ordinary refractive index of a droplet matches that of the scattering medium, 
then there will be no scattering of the polarisation component out of the plane of the incident 
wave and the uniform anisotropy in the droplet. In order to do this, we first extend algebraic 
relations that follow from Eqs. (46.1). These equations give the expansion of plane wave 
propagating in a uniformly anisotropic medium. We can rewrite them for the plane wave 
inside the droplet: 



A A A U AAA 

exp(ip e kj nc • f )£/ x (kj nc ) [Gj;(kj nc ) + - — sin 9i nc zj + exp(ip D kj nc • f )£/j / (kj nc ) ey(kj. 



j'm' 



n D ( exp(ip e k 

inc ' ^ e -^a; (kjn C ) Gy (kj nc ) exp(ip kj nc • T^jEy (kj^^Gj. (kj. 



(Al) 



(A2) 



where p e = m e p Q and m e = + u)/(l + u cos 2 6*j nc ) . The coefficients aj i " c ' ) and are 
defined by Eqs. (9.1) where the factor \ijn is changed to l/n . 

We can now combine the relations that come from definitions of the coefficient functions 
(see Eq. (45.1)) 



~ {inc) 



j'm 



j>>\m\ 



Y (a) *(f) 

jm \ I 



(A3) 



j'm' 



j'>|m| 



(a.e) / \ (inc) , (ct.m) / \ ~ (inc) 

n oqjy- m [Po) a),J + qj f . m J {po) ot),J 



r(a)*i 
jm 



(* ) " E [ 0°°> f ) + Qit (Po. f ) 1 ). , a G {m, e} , (A4) 



j'm' 



with the relations ( |A1|) - ([A2|) to evaluate the left hand side of the system (|3]) provided that 

r (c) ~ (c) •) r (inc) ~ (inc) 

i oi ■ ni ■ > — \ (y ■ ft ■ y 

V jmi "jmJ l"jm ) "jm J ■ 

To this end, we can use Eq. (|6]) to write down the sum on the left hand side of Eq. (|3|) in 
the following form: 
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o Pjj' 
(inc) (e,m) 
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i v^>"v/ \ ~(inc) 
j'm ' Qjj';m\P°) a j'm 
(inc) . (m,m) / \ ~ (inc) 



— n. 



2 if. ; m 

-l„( e >e) 

Fjj'-m 



(Po. 



j'm 
~ (inc) 

Or-, 

j'm 



(A5) 



It is seen that the elements of the column on the right hand side of this equation are the sums 
from the left hand sides of Eqs. (^3) and (|A4|) . On the other hand, from Eqs. (|A"l"l ) - (|A2|) , the 
square bracketed sums on the right hand sides of Eqs. (|A3])- (|A4j) are the plane waves. So, 
the elements of the column (|A5|) can be evaluated as scalar products of the vector spherical 
functions and the vector plane waves by using Eqs. (B6.I) and (B8.I) of appendix B in I. 
We can now apply this procedure to calculate the elements of the column (|A"5|) for the ordi- 



nary wave with {af™\ a/" c) } = {af™> y , 6cf™' y }- From Eqs. (p 
and E y (k inc ) = 1 in this case. The final result is 



H) we have E x (k inc ) = 



E 

j'>\m\ 



; m i 



(inc) 
^j'm; y 
~ (inc) 



/ Jj(po) \ 



a: 



(inc) 
jm;y 



n Q [jj(p )}' 





. -l ~(inc) 






n jj(po) 

\-\Jj(Po)]'J 



(A6) 



When n Q = n (and p = p a ), after multiplying ( |A~6|) by the matrices (|7|), we have 



((inc) \ / (inc) \ / (inc) \ /„\ 

- • E 1» " ffl • (A7) 



From these equations we immediately conclude that, when n = n and {«j™ c \ cuf^} 
{ a fwt> ^jnry}-i ^ ne solution of the system (H) is given by 



JLc) _ A inc ) a,( c ) _ ^( inc ) a( sca ) — o( sca ) _ n o( sca ) — ni sca ) _ n ( kR\ 

jm jm\y> jm jm;y> I jm rjm\y 5 / jm I jm;y ' V / 



So, the amplitudes of scattered wave flj™ y and pf^ y van i sri at n = n a . 
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